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We study a reaction-diffusion process that involves two species of atoms, immobile and diffusing. 
We start with only immobile atoms uniformly distributed throughout the entire space. Diffusing 
atoms are injected at the origin by a source which is turned on at time t = 0. When a diffusing 
atom collides with an immobile atom, the two atoms form an immobile stable molecule. The region 
occupied by molecules is asymptotically spherical with radius growing as t^^'^ in d > 2 dimensions. 
We investigate the survival probability that a diffusing atom has not become a part of a molecule 
during the time interval t after its injection. We show that asymptotically the survival probability (i) 
saturates in one dimension, (ii) vanishes algebraically with time in two dimensions (with exponent 
being a function of the dimensionless flux and determined as a zero of a confluent hypergeometric 
function), and (iii) exhibits a stretched exponential decay in three dimensions. 

PACS numbers: 05.40.-a, 82.20.-w, 81.65.-b 



I. INTRODUCTION 

Diffusion is the basic transport mechanism that un- 
derlies numerous phenomena in physics, chemistry, and 
biology JV^. In applications diffusing particles usually 
interact with each other, or with other species of parti- 
cles. Here we consider a system in which diffusing atoms 
(species A) are injected into a small localized region of 
a d— dimensional lattice. The entire lattice is initially 
occupied by immobile atoms (species one B atom 
per lattice site. Whenever an A atom hops to a lattice 
site occupied by a B atom, two atoms react to form an 
inert stable molecule (species i?*). Thus the process is 
described by the reaction scheme [3 |H] 

^(diffusing) -I- i? (substrate) i?* (stable) 

This model and its generalizations mimic important in- 
dustrial chemical processes such as electropolishing [7], 
dissolution of solids [5] , corrosion and etching [H IHl [H] j 
and erosion [10] . In the injection-controlled limit, namely 
in the situation when injection is so slow that each new 
A atom will undergo reaction before the next injection 
event, the process is also known as internal diffusion lim- 
ited aggregation (IDLA). 

The reaction A + B ^ B* proceeds at a certain rate 
which often greatly exceeds the hopping rate. Therefore 
we shall always assume that the reaction proceeds in- 
stantaneously, so each site contains either a B atom or 
a B* molecule. Hence as the process develops, the sys- 
tem can be separated onto two parts: The droplet that 
contains no B atoms (every lattice site inside the droplet 
is occupied by a B* molecule and can contain an arbi- 
trary number of A atoms) and the rest of the system that 
contains only B atoms. 

The droplet is of course a growing random set, but 
in the large time limit it becomes relatively more and 
more close to the ball. (This is not entirely obvious, e.g. 
external diffusion limited aggregates [11; are not spherical 
at all, their shape depend on the lattice and they have 
very intricate branching structures.) For IDLA, however. 



the assertion that the asymptotic shape of the droplet is 
a ball has been proved [T?; the general case when the 
strength of the source is finite has been also settled |13) . 
Intuitively, one expects [B] that in the long time limit the 
radius of the droplet R{t) and the density of A atoms 
inside the droplet, can be found by solving a Stefan-like 
problem, namely by solving a diffusion equation with a 
moving boundary whose position is determined in the 
process of solution [HI [15] . Gravner and Quastel proved 
[TFI that this is indeed correct. Using this reduction to 
the Stefan problem one finds [SJ [13] that the radius of 
the droplet scales as 



R oc 



^/t\nt d=l 
t^/'^ d>2 



(1) 



when t — > oo. 

Earlier simulations [51 [S] revealed that the surface of 
the droplets is strikingly smooth, e.g. in two dimensions 
the width of the interface appears to scale logarithmically 
with droplet radius. This logarithmic growth law has 
resisted the proof up until recently when it was shown 
(in the IDLA setting) that the width grows not faster 
than logarithmically [TBHT^ with the size of the droplet, 
more precisely not faster than In R when d = 2 [TTJ [18] 
and not faster than Vlni? when d > 3 [T71 [THj . 

In this article we examine first passage properties of 
diffusing atoms. One particularly interesting quantity is 
the survival probability S{to,t) that an A atom injected 
at a certain time to has not reached the (growing) outer 
boundary during the time interval (tg, t). We shall mostly 
focus on the limit when the observation time greatly ex- 
ceeds the release time: t ^ to. We shall show that in this 
limit the surviaval probability varies with the observation 
time t according to 



finite d = 1 

S (X It-'^ d = 2 

exp[— const X t^/'^] d = 3 



(2) 



where the exponent f3 depends on the (dimensionless) 
ratio of the flux rate to the hopping rate. 



2 



An outline of this paper is as follows. In Sect. [TTj 
we review the general formulation of the Stefan prob- 
lem and the solution of the Stefan problem in two di- 
mensions. The two-dimensional case is most common 
in applications as a localized injection from the 'third' 
dimension can naturally occur. The two-dimensional 
case is also mathematically more appealing as the Ste- 
fan problem admits a self-similar solution. In Sect. |III| 
we study the behavior of the survival probability in 
the two-dimensional setting. The analysis in the three- 
dimensional case is presented in Sect. |IV[ In Sect. [V| we 
returned to two dimensions, but assume that mobile A 
atoms are injected with a fluid, so their transport is de- 
termined both by convection and diffusion. Finally in 
Sect. [Vl| we give a brief summary. 



II. STEFAN PROBLEM IN TWO DIMENSIONS 

Consider a dimensional lattice initially occupied by 
immobile atoms (species B). Let diffusing atoms (species 
A) be injected at the origin. Whenever a mobile A atom 
lands on a lattice site with an immobile B atom, the two 
atoms immediately form a B* molecule: 



A + B ^ B* 



(3) 



The molecules B* are assumed to be immobile and sta- 
ble. The source of A atoms is turned on at time t ~ 0. 
The droplet that contains no B atoms is asymptotically 
a ball of radius R{t) and its boundary moves more and 
more deterministically when f S> 1. In the hydrodynamic 
framework, the concentration c(r, t) of A atoms (that is, 
the average number of A atoms per lattice site) satisfies 
a diffusion equation 



dc{r, t) 
dt 



1 d 



dr 



c{r,t) + J5(r)9{t) (4) 



inside the droplet < r < R{t). Here D is the hopping 
rate of A atoms, J is the strength of the source, and 9{t) 
is the Heaviside step function. 

The diffusion equation Q should be supplemented by 
the standard absorbing condition 



c(r = R{t),t> 0) = 
and the Stefan boundary condition 



dR 

lit 



-D 



, dc 
dr 



r=R 



(5) 



(6) 



which essentially expresses mass conservation (see [TUfT^ 
for derivations of such boundary conditions in various 
examples). The boundary moves and its position has 
to be determined in the process of solution. Therefore 
mathematically we arrive at the Stefan problem [Ml [15] . 

The original process occurs on the lattice and there- 
fore we set the lattice spacing to unity; this implies that 
the spatial coordinates r, the droplet radius i?(t), and 



the concentration c(r, t) are all dimensionless quantities, 
while the hopping rate D and the strength of the source 
have the dimension of inverse time: [D] = [ J] = 1/ (time) . 
Note that the ratio J / D is dimensionless and this param- 
eter plays an important role in the problem. It will prove 
convenient to use a slightly modified ratio 



J 

■inD 



which we shall call the dimensionless flux. 

We now focus on the two-dimensional case 
discuss the three-dimensional set-up in Sect 



IV 



(7) 



we shall 
In two 



dimensions, the Stefan problem (|4])-([6]) has been solved 
in [6 . Here we outline some details of the solution as we 
shall need them for the analysis of the survival probabil- 
ity. One seeks solution in the scaling form 



c(r,i) = c(o, e = 

The derivatives of the density read 



dc__R , 
dt R^"' 



dc 
dr 



1 , 

r' 



R 

Q^2 



i?2 



(8) 



(9) 



where (•)' = d{-)/d£^ and R = dR/dt. The governing 
diffusion equation Q reduces to 



c" + -c' 



RR 

IT 



(10) 



This equation is consistent if RR/ D, which is in principle 
a function of time, is actually a constant. Denoting this 
(yet unknown) constant by 2a we get 



R-" = AaDt 



and recast ( 10 1 into 

c" + (r'+20 c' = 



(11) 



(12) 



A solution to Eq. ( 12 ) satisfying the adsorbing boundary 
condition ([s]) reads 



m=Af'^e- 



(13) 



Plugging this solution to the Stefan boundary condition 
^ allows one to express the amplitude A via a: 



A = 2ae" 
To determine a, we notice that 



Jt= I c{r, t) 27rr dr + ttR^ 
Jo 



(14) 



(15) 



Equation ( 15 ) is just the conservation of the total number 



of A atoms. Re- writing ( 15 1 through the scaling variables 

(16) 



and using (111 we arrive at relation 



a e 
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which implicitly determines a in terms of the dimension- 
less flux $. Note that the fraction of diffusing A atoms 
(that is, those which are not the part of i3* molecules) 
with respect to all emitted A atoms is 1 — e~". Put differ- 
ently, the total number A(t) of diffusing A atoms grows 
linearly with time: A(t) ~ (1 — e^°')Jt. 



III. SURVIVAL PROBABILITY IN TWO 
DIMENSIONS 

Although the Stefan problem described in Sect. [lT| has 
been previously understood as far as the general growth 
characteristics are concerned, more subtle features char- 
acterizing the behavior of diffusing atoms haven't been 
explored. Here we analyze first passage properties [50] 
of diffusing atoms. As an example, consider an A atom 
that has been injected at a certain time to a-nd ask for the 
survival probability S(to,t), namely the probability that 
it has not reached the (growing) outer boundary during 
the time interval {to,t). The two- variable survival prob- 
ability S{to, t) apparently has an interesting scaling form 
in the limit when both and t are very large and com- 
parable, that is, 



in 



t 



finite 



(17) 



More simple behavior should arise when the observa- 
tion time greatly exceeds the release time, that is, t ^ to. 
Indeed, the average distance from the origin exhibits the 
diffusive growth, (r^) — AD{t — to), which is similar to 
the growth of the radius of the droplet, R = y/ADat. 
This suggests that the survival probability exhibits an 
algebraic decay 



S{tQ, t) (X t 



-13 



when t ^ to 



(18) 



The persistence exponent /3 should be a function of a. 
The dependence /3 = I3{a) can indeed be analytically de- 
termined. We shall employ essentially the same method 
as in Refs. [KMill]- 

To determine the survival probability (in the t ^ to 
limit) we look at a more detailed quantity, the proba- 
bility distribution P{r,t). By definition, P{r,t)dr is the 
probability that the atom has never reached the growing 
outer boundary of the droplet during the time interval 
(to,t) and at the final time the atom lies in the region 
of area dr around r. One anticipates that the probabil- 
ity distribution P{r, t) is actually a function of the scaled 
spatial coordinate ^ = r/R{t). Therefore we make a scal- 
ing ansatz similar to (l8| 



m 



(19) 



We shortly write S{t) instead of S{to,t) since we assume 
that to ^ t and we are mostly interested in the depen- 



dence on t. The time-dependent pre-factor in ( 19 1 has 



been chosen to ensure the validity of the connection be- 
tween the probability distribution P(r, t) and the survival 
probability. Indeed, 

S{t)^ [ P{r,t)2nrdr = S{t) [ dC^T(0 
Jq Jq 

and therefore the consistency condition is 



(20) 



The probability distribution P(7', t) satisfies the diffusion 
equation 



dP 
~dt 



= D 



1 d 



(21) 



in the region < r < R{t). There is no source, apart from 
the fact that initially the atom is at the origin. Hence 
the initial condition is 



P{r,to)^S{r) 
and the adsorbing boundary condition is 
P{r = R{t'),t') = 



(22) 



(23) 



for to < t' < t. In the t 3> to limit, the survival proba- 
bility exhibits a power-law dependence ^(to, t) cx t~^ on 
the observation time, and as long as we are interested in 
the persistence exponent /3 we can forget on to. 

Using (19) and (11) we recast the diffusion equation 



(21 1 into an ordinary differential equation 

y" + (r^ + 2aO y + 4a(l + /3)J' = 
The change of the variable ^ to 

c=-ae 



transforms ( 24 ) into a hypergeometric equation 

dy 



d^y 



{i + l3)y = Q 



(24) 
(25) 

(26) 



The solution to ( 26 ) must satisfy the adsorbing boundary 



condition ( 23 1 , or equivalently 

y(e = 1) = 0, 



(27) 



and the integral requirement (20). It should also be reg- 
ular at r = 0. 

The condition of regularity at r = selects a solution 
up to an amplitude, y — C F{1 + f3, 1; C), where F is the 
confluent hypergeometric function |23j . Therefore 



y ^ c F{i + 



(28) 



and the boundary condition ( 27 ) yields an equation 
F(H-/3,l;-a) = 



(29) 
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which relates the persistence exponent /? and the growth 
constant a. 



The /3 — /3(a) dependence provided by Eq. (291 is im- 
phcit. This equation has many solutions, the proper one 
corresponds to minimal /?. Having determined (3 — /3(a), 
we can fix the amplitude C in Eq. ( 28 ) by requiring the 
validity of ( 20 ) . The final result reads 



y = 2 



Fil + p,l;-ae) 
F(l + /3,2;-a) 



(30) 



More explicit results can be established in certain spe- 
cial cases. For instance, since /3 — >■ oo as a — >■ 0, we 
can asymptotically re-write the confluent hypergeomet- 
ric function F = F{1 + /3, 1; —a) in terms of the Bessel 
function of zero order: 



S{to,t) on the release time <o- On dimensional grounds, 
one anticipates a simple algebraic behavior 



Sito,t) 



(34) 



Equation ( 34 ) is expected to hold in the same t ^ limit 



as (18); in contrast to the latter, however, Eq. (34) has a 



proper dimensional form. Moreover, one anticipates that 
the survival probability remains finite in the double scal- 
ing limit (17) and Eq. (34) agrees with this requirement. 



This suggests that Eq. (34 1 gives an approximately cor- 



rect behavior in the entire time range. 



IV. THREE DIMENSIONS 



F = 



1 - 
1 - 



(l + /3)a (l + /3)(2 + /3)a^ 

(1!)2 
/3a (/3a)2 



(2!) 



IA2 



(1!)2 (2!)2 
= Jo(2y^) 



(31) 



The Bessel function has infinitely many zeros, Jo{x) — 
when X — ±xi,±X2, ■ ■ ■■ The relevant root is closest 
to the origin, namely xi = 2.404825558.... Therefore 
/3 = (xi/2)Va = 1.445796491/a in the a ^- limit. 
When a is small, the dimensionless flux is also small, 
<f> ~ a according to Eq. ( 16 ) , and hence 



1 445796491 D 
^ L.^^!,iymyi ^ 18.16841454 - 

$ J 



(32) 



when D ^ J. In this limit the probability distribution 
also acquires a neat form 



Ji{xi) 



(33) 



The confluent hypergeometric function simplifies when 
the indexes are integer. For instance, let /3 — 1. Since 
F(2,l;-a) = (1 - a)e-", equation ([29]) yields a = 1. 



Recalling ( 16 ) we conclude that 



when $ = e 



and the probability density in this case is 
a' = 2(l-e^)e^~^' 



For ^ = 2, we have F(3, 1; -a) = {l ~ 2a + la'^)e- 



Relation (|29|) leads to 1 — 2a -f ^a — 0, from which 



a = 2 — y/2 and therefore 



when $ 



/2) 



Previous scaling analysis does not allow one to probe 
the dependence of the two-variable survival probability 



We begin again with the Stefan problem. In three di- 
mensions, the governing diffusion equation reads 



dc 

di 



= D 



2 _9 
r dr 



(35) 



and we seek a solution to ( 35 1 inside the growing droplet 



< r < R{t). The boundary conditions are the same as 
in two dimensions, E qs. (|5| )-(|6|. Using again the scaling 
ansatz ([s]) we reduce (|35| to 



2 , RR , 

— c = fc 

£, D ^ 



The conservation law 
Jt 



c{r, t) Anr"^ dr 



47r 



(36) 



(37) 







tells us that Jt> ^R^ 



Hence the radius grows at most 
as t^^"^ and RR decays at least as i"^/'^. Therefore the 



right-hand side of equation ( 36 1 asymptotically vanishes 
and we get c" -|- | c' = from which 



(38) 



Near the origin the density is stationary and it satisfies 
V^c = — -^(5(r). Therefore c — ^/r when r <^ R. Com- 
paring this asymptotic with (38) wc find C = ^/R. This 
allows us to re- write ( 38 ) as 



c = $ 



1 



Note that the number of diffusing A atoms is 

ni? 



A{t) 



c(r, t) Airr dr 







J i?2 

D 6 



(39) 



(40) 



The number of bounded A atoms (which are included 
into B* molecules) scales as R^ and therefore it consti- 
tutes the dominant part of all injected A atoms. Hence 



5 



Jt 



-B? implying that asymptotically the radius of 



the droplet is given by 



R = 



47r 



1/3 



(41) 



Note that the fraction of injected A atoms which continue 
to diffuse, (Jt)-U(i), decays as i(3$)2/3(7:»t)-i/3. 

Since the radius grows slower than diffusively, the sur- 
vival probability should decay very quickly. To establish 
this decay qualitatively we must solve the diffusion equa- 
tion 



dP 
~dt 



d(^ 



2 d_ 
r dr 



P 



(42) 



subject to the initial-boundary conditions ( 22 1-( 23 ). The 



behavior is again simpler than in two dimensions. The 
leading asymptotic is 



P{r,t)^f{t) 



sin(7rr/i?) 



(43) 



where sin(7rr/i?) is the eigenfunction of the Laplace 
operator corresponding to the smallest eigenvalue. The 
dominance of the eigenfunction corresponding to the 
'ground state' is_easy to appreciate — since the growth 
is slow, R <^ V Dt, the density has enough time to equili- 
brate, so its spatial behavior is the same as in the case of 
the fixed-size droplet and only the time-dependent fac- 
tor / = f{t) may differ. For more details (in the one- 
dimensional setting) see [21]. 
Using ( 43 1 we find 



1 dP 
P Ik 



/ 
/ 



1 



V'P 



TTT R 

^ R 

TT 

R 



(44a) 
(44b) 



It turns out that the first term on the right-hand side of 
(44a I is dominant. Hence f = —D(j^) f from which 



f{t) - exp<^ -D / dt' 



R{t') 



(45) 



Equation ( 45 ) is valid for arbitrary growth law of the 
radius which is slower than diffusive, R <^ %/ Dt, e.g. for 
any algebraic growth R ^ t°- with a < 1/2. Specializing 



the general result ( 45 ) to our case when the radius grows 
according to Eq. (41 ) we obtain 



/~exp{-^2^-2/3 (3^i)i/3| 



(46) 



Note that /// - t~2/3^ ^Yiile R/R - . Hence the first 
term on the right-hand side of (44a) is indeed dominant. 



Finally, the survival probability is given by 

„2 , 



S = 



P{r,t) Attt dr 



4/i?2 



du u sin u 



nr/R) 



The term 4i?^ provides just a power-law correction to a 
controlling exponential decay of /. We have ignored such 
terms in a derivation of (45 )-(46 ). Therefore the survival 



probability decays according to the same law (46), viz. 



S ~ exp-^ 



-2/3 



(47) 



confirming the announced result ^ in three dimensions. 

In the three-dimensional setting, one can analytically 
determine the asymptotic behavior of the two-variable 
survival probability Sitg, t) even when the realize time 
is not negligibly small in comparison with the observation 
time t. The probability density P{Q,to;r,t) describing 
this situation can be expressed using a quasi-static ansatz 
similar to ( 43 1 , namely 



P(0,to;r,t) = /(io,t) 



sin(7rr/i?) 



(48) 



This description is applicable when the typical diffu- 
sion length characterizing the released atom, a/ D{t — to), 
greatly exceeds the droplet radius R{t). Therefore using 
( 41 1 we arrive at the criterion of the validity of ( 48 1 



2/3 



(49) 



Using ( 48 1 instead of ( 43 ) we find the controlling expo- 



nential behavior of the amplitude in Eq. ( 48 ) 
fito,t) ~ expJ ~D dt' 



R{t') 



(50) 



which gives again the controlling exponential behavior of 
the two- variable survival probability 



S{to,t) exp 



$2/3 



{3Dty/^ - {3Dto 



,1/3 



(51) 



It is tempting to use Eq. (51 1 to compute the number 
of diffusion A atoms 



A{t) = J dtoSito,t) 



(52) 



and to compare the outcome with the asymptotically 
exact prediction (40). The asymptotic results match if 
we use (51 1 with numerical factor ^tt^. This seemingly 



gives us the exact numerical amplitude and, more impor- 
tantly, tells us that there is no power-law pre-factor to 
the prediction (51). A closer examination shows, how- 



ever, that the main contribution to the integral on the 
right-hand side of Eq. (52 1 is gathered in the region 
-i/3(i/(j))2/3^ which is precisely the temporal 



D- 



t-to 

range [see (49 )] where the starting ansatz (48 ) is no longer 
valid. Therefore we cannot use the sum rule (52) to fix 



the correcting pre-factors to the controlling exponential 
behavior, although it does seem plausible the power-law 
pre-factor is absent. 
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V. REACTION-CONVECTION-DIFFUSION 
PROCESS WITH A LOCALIZED SOURCE 

In this section we investigate what happens if diffusion 
is supplemented by source-driven convection. The flux of 
diffusing A atoms can be organized through the flux of 
the fluid containing A atoms. Here we consider the most 
natural two-dimensional situation. We thus assume that 
there is a flux of both A atoms and fluid; as before, we 
posit that B's and B*^s are attached to the lattice and 
therefore immobile. This convection-diffusion problem 
may be interpreted as a pure diffusion problem [24j in 
a space with an effective dimension different from the 
physical dimension dphys = 2. More precisely, let Q be 
the fluid flux; put differently, the velocity field is given 
by y = 7p-. It turns out that this convection-diffusion 
problem can be recast as a pure diffusion that occurs in 
a space of effective dimension [24] 



d = 2-2q, q = 



Q 

4ttD 



(53) 



Therefore the analysis is similar to the one described in 
the previous sections. Most of the results which we will 
derive below hold for both sink and source flows, but 
the case of source flows {Q > 0) is more natural, so we 
shall assume that we have a source of both A atoms and 
fluid. Note that for source flows g > 0, so the effective 
dimension satisfies d < 2. 



A. Stefan Problem 

Generally in two spatial dimensions the convection- 
diffusion equation for radial velocity field V = V{r) reads 



dc 
di 



V 



dc 
dr 



D 



1 d_ 
r dr 



(54) 



For the velocity field V 
equation becomes 



— 7^ this convection-diffusion 



dc 

at 



D 



Qj.2 



1 



2q d_ 
dr 



c + JS{r)9{t) (55) 



where we have added the source term. The comparison of 
( [SSj ) with Q explains the assertion ( 53 ) that the effective 
dimension is equal to d = 2 — 2q. The radial symmetry 
and the form of the velocity field are of course es- 
sential for the above reduction of convection-diffusion to 
pure diffusion j25j . 

Using the procedure detailed in Sect. Ill] we obtain 



■ l-2g 



2a^) c' = 



instead of (12 1, from which 



c(0 =^ / drjr] 



(56) 



(57) 



The amplitude A and the constant a are related through 



the same relation (14 1. Using again the conservation law 



(151 in conjunction with (II) and (57 1 we obtain 



a + Aa dr]r]^+^ie-°"^^ = $ 

"'0 



(58) 



When q = we recover our previous results. Here are 
two more examples. 



1. Effectively zero- dimensional case 

When q = 1, or equivalently d — 0, the integral on the 
left-hand side of equation ( 58 1 is elementary and ( 58 ) 
simplifies to 



e" - I = $ 



(59) 



This effectively zero-dimensional case is easier than the 
original two-dimensional case as the density (57) also sig- 
nificantly simplifies: 



c(0 = (1 + <f ) 



i-e 



- 1 



(60) 



2. Effective d = —2 dimension 



When q = 2, or equivalently rf = — 2, equation (58) 
simplifies to 



„ e" - I - a ^ 

2 = $ 

a 



(61) 



while the density (57) becomes 



c(0=e'e"(i-«-)-l + ^ (62) 

a 



B. Survival Probability 



Assuming the validity of the scaling ansatz ( 19 1 we 
obtain an ordinary differential equation 

T" + (i^+2ae)T'+4a(l-f/3)T = (63) 



for the scaled probability distribution. Transforming 
again the variable C to C which is defined by Eq. (251, 
we recast ( 63 ) into 



d^T dT 

C^ + (l-9-C)^-(l + /?)y-0 (64) 

This hypergeometric equation has two linearly indepen- 
dent solutions, 

F(l + /3,l-9;-<') & e^Fil + P + q^l + q'-ae) 
and for g > the latter solution is selected: 

y = C e'«i^(l + /3 + g, 1 + g; -<') (65) 
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The boundary condition (27) yields equation 



F{l + l3 + q,l + q:-a)=0 



(66) 



The normalization requirement ( 20 ) allows one to fix the 



amplitude C in ( 65 1 to yield 



T-(2 + 2«)£2'^^^ii±^±^ll±^i^ (67) 
■'~^' + "^>^ F(l + /3 + g,2 + g;-a) ^''^ 

More explicit results can be obtained in the small flux 
limit. In this situation a ~ $ as it follows from (58 1, and 
when a the persistence exponent diverges /? — > oo. 
In the limit 



a 



0, f3 



oo, 



finite 



the same computation as in ( 31 1 allows one to express the 



confluent hypergeometric function via the Bessel function 
F{f3,l + q;~a)o,r{q + l):^^l^y^ (68) 



Comparing (66) and (68) we get 2 -/a/3 — xi^q, where 
xi^q is a root of the Bessel function of order q, that is 
Jq{xi,q) = 0; there arc infinitely many such roots and 
Xi^q is actually the smallest positive root. Thus in the 
a — > the persistence exponent diverges as 



(69) 



Using the asymptotic form ( 68 ) we find that the proba- 



bility distribution ( 67 1 simplifies to 

Xl,q ^''Jq{xi,q£,) 
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Jl+q{xi^q) 



(70) 



The physical requirement that the probability distribu- 
tion must be non-negative explains why the proper root 
Xi^q is the smallest positive root of the Bessel function. 



The probability distribution (67) also simplifies when 



/3 is integer. For instance, when P — 1 we get a — 1 + q. 
In this case 

s ^ t-\ T = 2(1 + q)e'{i - e) e(i+')(i-«') 

As an example, consider effectively zero-dimensional 



(661 become 



case. When g = 2, or equivalently d — Eqs. (65) and 
T = Ce'^^(2 + /3,2;-<') 



and 



F{2 + l3,2--a) = 



(71) 



(72) 



When a — )■ 0, the persistence exponent diverges /3 — )• oo. 
The same computation as in (31 ) shows that in this limit 



which determines the persistence exponent (3 — f3(q,a). order 



the confluent hypergeometric function that appears in 
( 72 1 can be expressed via the Bessel function of the first 



F(2 + /3,2;-a) 



Therefore in the a — ^ limit 



P = ^^i^ = 46.12477111 ^ 
a J 



(73) 



where xi — 3.83170597 is the smallest positive root of 
the Bessel function of the first order and in deriving the 
second equality we have used the asymptotic $ ~ a. In 
this limit the probability distribution ( 71 ) also simplifies 
and acquires a neat form 



xiS,Ji{xii) 

J2{xi) 



(74) 



VI. SUMMARY 

We have studied the behavior of the two-variable sur- 
vival probability 5(^0, t)- We have mostly focused on the 
dependence on the observation time t in the limit when it 
greatly exceeds the release time, t ^ tg. We have demon- 
strated that in two dimensions the survival probability 
exhibits an algebraic decay (18 1. We have also shown 



that the persistence exponent is related to the dimen- 
sionless flux via a root of the confluent hypergeometric 
function. Similar results continue to hold when mobile 
A atoms are injected with a fluid, so their transport is 
determined both by convection and diffusion. 

In three dimensions, we have shown that the survival 
probability exhibits a stretched exponential decay (47|. 
Furthermore, we have computed the asymptotic behav- 
ior of the two-variable survival probability for (almost) 
arbitrary release and observation times. 

We haven't yet derived the announced result ^ in one 
dimension. To fill this gap we recall the well-known result 
(see e.g. Refs. (SI USD that in one dimension the droplet 
grows a little bit faster than diffusively, namely there is 
a logarithmic correction R = [2Dt ln( J^t/D)]^/^ to the 
diffusive growth. Using e.g. the approach of Ref. [3T] 
one finds that the survival probability remains finite in 
the t — oo limit; in other words, there is a finite chance 
that an A atom will never meet a B atom, that is, it will 
forever remain inside the growing droplet. 
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